Perturbation Theory in Lemaitre-Tolman-Bondi Cosmology 
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The Lemaitre-Tolman-Bondi solution has received much attention as a possible alternative to 
Dark Energy, as it is able to account for the apparent acceleration of the Universe without any 
exotic matter content. However, in order to make rigorous comparisons between these models and 
cosmological observations, such as the integrated Sachs- Wolfe effect, baryon acoustic oscillations 
and the observed matter power spectrum, it is absolutely necessary to have a proper understanding 
of the linear perturbation theory about them. Here we present this theory in a fully general, 
and gauge-invariant form. It is shown that scalar, vector and tensor perturbations interact, and 
that the natural gauge invariant variables in Lemaitre-Tolman-Bondi cosmology do not correspond 
straightforwardly to the usual Bardeen variables, in the limit of spatial homogeneity. We therefore 
construct new variables that reduce to pure scalar, vector and tensor modes in this limit. 
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I. INTRODUCTION 

The concordance model of cosmology is a result of as- 
suming that the Universe is approximately the same at 
all points in space, and that at every point it appears 
isotropic. The observations we make locally can then be 
extrapolated to the Universe at large, and the resulting 
Friedmann-Lemaitre-Robertson- Walker (FLRW) model 
can be used to interpret all cosmological observations. 
The FLRW models are, of course, very well understood, 
and have proved to be highly successful in very many 
different capacities. Unfortunately, despite its aesthetic 
appeal and mathematical simplicity, the adoption of the 
FLRW cosmology has led to the undesirable requirement 
that the Universe must be filled with a smoothly dis- 
tributed and gravitationally repulsive exotic substance, 
known as Dark Energy. The conceptual problems asso- 
ciated with Dark Energy are well known, and are dis- 
tressing enough that it now seems worth investigating 
alternative, less appealing models. 

One such alternative is the spherically symmetric, dust 
dominated cosmological solution of Lcmaitrc, Tolman 
and Bondi This model is isotropic about a single 
point in space, in agreement with local galaxy surveys 
and observations of the Cosmic Microwave Background, 
but is anisotropic everywhere else 0. As such, we are 
required to live very near the centre of symmetry, indi- 
cating a break with the Copernican Principle that "the 
Earth is not in a central, specially favoured position" 
Q. The Lemaitre-Tolman-Bondi (LTB) model allows the 
Hubble rate to vary with radial distance from the cen- 
tre of symmetry, as well as with time, and so allows for 
the possibility of apparent acceleration (if observations 
within them are interpreted in an FLRW frame- work), 
while no part of them need undergo any accelerating ex- 
pansion. 

One may then try and argue on aesthetic grounds, or 
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on the grounds of presumed early Universe physics, which 
of these models is more likely. Clearly there are good rea- 
sons to disfavour either, depending on which one abhors 
most: the Earth being in a special position in the Uni- 
verse, or the existence of Dark Energy, which requires 
us to live at a special time when the acceleration is just 
starting to kick in. However, a more objective (and sci- 
entific) approach to the problem would be to try and 
distinguish between the two models with direct observa- 
tions. This has been attempted a number of different 
times in the literature already Q. However, in order to 
make comparisons to different types of astrophysical data 
it is absolutely necessary to have a proper understand- 
ing of how structures form and behave in these models, 
and this means understanding their perturbation theory. 
An alternative approach is to try to test the Copernican 
principle directly [5j, but this is very difficult to achieve 
at present through direct tests. 

For observations such as the integrated Sachs- Wolfe ef- 
fect, Rees-Sciama effect, or matter power spectrum it is 
obvious that a proper interpretation of the data is depen- 
dent on understanding the relevant perturbation theory. 
However, other observations, that do not necessarily re- 
quire perturbation theory in FLRW, will also require a 
knowledge of perturbation theory in the LTB case. An 
example is the use of baryon acoustic oscillations as a 
probe of cosmology. In FLRW it suffices to consider co- 
moving points in an unperturbed fluid, to see how length 
scales at recombination are evolved onto distance mea- 
sures on our past light cone. In an LTB universe, how- 
ever, the rate of growth of structure can be different at 
different places, and the peculiar velocity of perturba- 
tions need not cancel when averaged over suitable scales 
(in fact, one may naively expect otherwise). Likewise, 
when considering light reflecting off a galaxy cluster in 
FLRW it may suffice to consider that on average the mo- 
tions of galaxy clusters will cancel out, so that approxi- 
mating galaxy clusters as being comoving with the back- 
ground fluid may be sufficient to get accurate results. In 
an LTB universe this is unlikely to be true, as gradients 
in gravitational potentials will develop due to different 
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growth rates at different positions. 

As well as the requirement of understanding linear per- 
turbation theory in order to accurately account for al- 
ready established probes of cosmology, there is also the 
possibility of finding new observational phenomena that 
could be used to distinguish the LTB and FLRW cos- 
mologies, or constrain the type of LTB models that are 
observationally viable. An example of this is the poten- 
tial coupling of scalar gravitational potentials to tensor 
gravitational waves. In an FLRW background (to linear 
order) the decomposition theorem ensures that no such 
coupling occurs. On a homogeneous background scalar, 
vector and tensor modes then evolve independently. In an 
inhomogeneous LTB background, however, the decompo- 
sition theorem will only be valid on the two dimensional 
surfaces of spherical symmetry. In the three dimensional 
space there will be no such theorem, and in general one 
should expect coupling to occur between what we conven- 
tionally think of as 'scalar', 'vector' and 'tensor' modes. 
Though complicated, such effects could potentially be 
used as probes of the geometry of the universe. 

In section [TT] we introduce the formalism required for 
perturbation theory about a spherically symmetric back- 
ground. In section UlI Al and IIII Bl we present the pertur- 
bation equations, and in section UlI CI and IIVI we discuss 
their solutions and how they relate to perturbation the- 
ory in FLRW. Finally, in section [V] we conclude. The 
Appendices contain information about gauge invariant 
quantities for the metric perturbation, and the perturba- 
tion equations for some special cases. 

We use units in which c = G = 1, throughout. 



II. PERTURBATION THEORY ON A 
SPHERICALLY SYMMETRIC BACKGROUND 

Perturbations on a spherically symmetric space-time 
have been considered a number of times in the literature, 
mostly in the context of modelling static and stationary 
stars [6] . A time dependent formalism was first developed 
by Seidel Q, and later specialised to a space-time con- 
taining a perfect fluid by Gundlach and Martm-Garcia 
(GMG) Perturbations of self similar models have 
been investigated in Q, and have the advantage that 
all PDEs reduce to ODEs in the system of perturbation 
equations. A limited class of perturbations in LTB cos- 
mology have also been studied by Zibin [10( . 

Different formalisms exist in the literature for con- 
structing both the system of perturbation equations, and 
the gauge-invariant variables. Gerlach and Sengupta [llj 
(GS) developed a formalism based on a 2+2 covariant 
split of the metric and energy-momentum tensor. This 
turns the field equations into a system of second-order 
PDEs, and was made explicit by GMG, who found a 
closed set of master equations describing all perturba- 
tions. Their approach was used in the study of self- 
similar LTB perturbations performed in Q. Alterna- 
tively, a covariant 1+1+2 formalism has been developed 



in [12j, which builds on the covariant 1+3 formalism that 
has been usefully applied in cosmology [Hj]. Here the 
Bianchi and Ricci identities, plus the Ricci rotation co- 
efficients for the semi-tetrad introduced, are covariantly 
split into a system of first-order differential equations. 
This formalism has not yet been reduced to a tractable 
set of master equations for a general space-time, but has 
been used by [lCj] to study LTB perturbations. Here we 
apply the GS formalism developed by GMG to the case 
of an LTB space-time. This results in a simple set of 
coupled second-order PDEs that describe general pertur- 
bations to the space-time, in a gauge invariant way. 



A. The LTB background 



The unperturbed LTB line-element can be written 



d S ^ 



a\{t,r) 
(1 — nr 2 ) 



dr z +a 2 L {t,r)r z dtt\ (1) 



where a\\ = (ra±) tr and k = n(r) is a free function of 
r. The FLRW scale factor, a, has been replaced here 
by two new scale factors, ay and a±, describing expan- 
sion parallel and perpendicular to the radial direction, 
respectively. For future use we will define the radial and 
azimuthal Hubble rates to be 



an aj_ 
H II = — and H± = — , 

an a 



(2) 



where an over-dot denotes partial differentiation with re- 
spect to t. The analogue of the Friedmann equation in 
this space-time is then given by 



2 M K 
Hj_ = — s o-, 



(3) 



where M = M(r) is another free function of r, and the 
locally measured energy density is 



87rp 



(Afr 3 ), r 



which obeys the conservation equation 
p+(2H ± + Hu)p = 0. 



(4) 



(5) 



The acceleration equations in the perpendicular and par- 
allel directions are 



ax M , ail M 

— = — — — o~ and — = — 47T/9 H — s-. (6) 

a± 2a] an a] 



For what follows it will also be useful to define the radial 
derivative 



(7) 
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This derivative does not commute with the time deriva- 
tive, but instead obeys 



{X)'-{X')=H l X'. 
We also define the curvature function 



W EE 



a±r 



The following relations are then obeyed 
H' ± = W(H {] -H x ), 
W = -H X W, 

W = —W 2 - Airp + HxHn - 



M 
2^' 



(8) 
(9) 

(10) 
(11) 
(12) 



In the perturbation equations that follow we will choose 
to eliminate k in favour of W, so that the equations take 
their simplest form. 



B. Harmonic functions 

A natural way to split this space-time is in a 2+2 de- 
composition, so that the space-time manifold becomes 
M 4 = M 2 x S 2 , where S 2 indicates the 2 dimensional 
spherically symmetric surfaces. We will use lower case 
Latin indices a, b, c, ... to denote coordinates in S 2 , up- 
per case Latin indices A, B, C, . . . to denote coordinates 
in M 2 , and Greek indices ^, v, £, . . . to denote coordinates 
that run over all 4 space-time dimensions. 

In FLRW cosmology, any perturbation can be split 
into scalar, vector and tensor (SVT) modes that decouple 
from each other, and so evolve independently (to first or- 
der). This classification is based on how they transform 
on the homogeneous and isotropic spatial hyper-surfaces, 
and is essentially just a generalisation of Helmholtz's the- 
orem [Til ]. 

Such a split cannot usefully be performed in the same 
way here, as the background is no longer spatially ho- 
mogeneous, and modes written in this way would couple 
together (as we shall see). However, one can perform an 
analogous classification based on how the perturbations 
transform on the surfaces of spherical symmetry. This 
results in a decoupling of the perturbations into two in- 
dependent modes, called 'polar' (or even) and 'axial' (or 
odd), which are analogous, but not equivalent, to scalar 
and vector modes in FLRW. Unlike the FLRW case, how- 
ever, there is no further decomposition into tensor modes 
as no non-trivial symmetric, transverse and trace- free 
rank 2 tensors can exist on S 2 . Therefore, only two dis- 
tinct sectors exist. Scalars (rank tensors) on S 2 can 
then be expanded as a sum of polar modes, and higher 
rank tensors on S 2 can be expanded in sums over both 
the polar and axial modes. Only a scalar can contain 
spherical perturbations (given by £ = 0, defined below), 
and only scalars and vectors (i.e., tensors of rank and 
1) can contain a dipole term (£ — 1). Higher multipoles 
can be present in all tensors. 



An appropriate family of basis functions for this split 
are tensor spherical harmonics. These are derived from 
the usual spherical harmonic functions, Y^ m \x a ), that 
obey 



(13) 



where the Natural number, £, gives the angular scale of 
the perturbation. The Laplacian, V 2 here, is on the 
surface of spherical symmetry, and is given by V 2 = 
7 ab ^:afc, where the colon subscript indicates a covariant 
derivative with respect to the metric on the unit sphere, 
jab- Scalar perturbations on S 2 can then be written with 
their angular dependence given in terms of the solutions 
to this equation, and expanded as 



oo rn—t 
=0 m=-l 



tm \x A )Y {tm \x a ). (14) 



It is now possible to construct a basis for all higher 
rank tensors from Y^ m \ its covariant derivatives, and 
the contractions of those derivatives with the fundamen- 
tal antisymmetric tensor, e ao . Modes that can be de- 
scribed without requiring e a b are called polar, while those 
that require e a b are called axial. 

We can now form harmonic functions for higher rank 
tensor perturbations with polar degrees of freedom by 
first defining the vector, for £ > 1: 

yj' m ) = Y. { a tm \ (15) 

We define also the trace-less tensor, for £ > 2: 



Y 



Y\ 



+ 



'lab- 



(16) 



Taking divergences of Y^ 1 ™^ and Yj£ m ' reduces these ex- 
pressions to equations involving y( £m ). 

For axial perturbations on S 2 we define a divergence- 
free vector harmonic, for I > 1, 



b v {tm) 
a 1 :b 



(17) 



We can then construct a symmetric and trace-free rank-2 
axial harmonic function by defining, £ > 2, 



Y, 



(a:b) 



= -2e a (n Y. 



(JLr, 



(a I :b)d ' 



(18) 



where round brackets around indices denote symmetri- 
sation, as usual. Taking covariant derivatives, this ten- 
sor harmonic can be reduced to an expression involving 

1 a 

For both parities, the 'vector harmonics' obey 
(similarly for Y„ f m ^), and the 'tensor harmonics' obey 



v 2 r> 



[4-£(£+l)]Y^ 



(20) 
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Because the vector harmonics are orthogonal for each 
£, any rank-1 tensor perturbation can now be expanded 
as 



b a (x A ,x a )=J2 £ ^ im \x A )Y^ m \x a ) 

4>^ m \x A )Y^ m \x a ), (21) 



t=l m=- 



where <^ im > and 0^ m ) are independent for each (£m), 
and are given by 



lm) _ 



(22) 



and 



— J dn(e a ^ b .. a )Y { * em) . (23) 



£(£■ 



There are no I = vector dof because £ = describes 
spherical modes. 

Finally, let us consider trace-less rank 2 tensor pertur- 
bations on S 2 . A suitable orthogonal basis in this case, 
for £ > 2, is given by the two harmonic functions Y^v 

and Y^ m \ defined above. Any rank 2 trace-free tensor 
perturbation can then be written 



oc m=l 



1=1 m =- 



l V)^"V)> (24) 



where 



V (£ + 2)\ 

M- 2)' 



b 



'' !: " , ' 2 ^T^J An ^a C r.. bc )Y ( \ m) . (25) 



As our background is spherically symmetric, all per- 
turbations with different (£m) decouple from each other. 
This is analogous to the decoupling of Fourier modes in 
FLRW. We shall therefore drop the (lm) labels on all 
quantities that follow. Where a function of x A is mul- 
tiplied by a harmonic function, a sum over £ and m is 
implied. 



C. Metric perturbations 

Gauge invariant variables for general perturbations of 
a spherically symmetric background, with arbitrary mat- 
ter content, have been formulated by Gerlach and Sen- 
gupta (GS) [III- We reiterate their results, relevant to 
the present study, in Appendix [A] We then review the 
perfect fluid formalism of GMG, specialised to a dust 
filled universe, in Appendix [Bj These studies show that 
there is a preferred gauge, known as the Regge- Wheeler 



(RW) gauge [15| , in which the perturbation variables are 
equal to gauge invariant quantities (comparable to the 
longitudinal, or conformal Newtonian gauge, of FLRW). 
The RW gauge is the choice that any off-diagonal polar 
modes in the metric that have an angular index are zero, 
and that in the axial modes there is no perturbation to 
the g ab components (i.e. the S 2 components). Calcu- 
lations performed in the RW gauge are then equivalent 
to those performed when considering the gauge invariant 
variables of Q and [Tlj |. 

The general form of polar perturbations to the metric 
can now be written, in RW gauge, as 1 

ds 2 = - [1 + (2ry - x - ¥>)Y] dt 2 - o dtdr (26) 



a 2 dr 2 



+ [1 + (X + <p)Y] —J + a 2 ± r 2 (l + ^Y)dtf , 

(1 — KT J 

where ry(i, r), x(i,r), ip(t,r) and <j(t, r) are equal to the 
gauge invariant quantities of GS and GMG, as shown in 
Appendices and [B] The general form of polar matter 
perturbations in this gauge are parameterised by 



P 



Ua + \ "'I! A 

LTB i 



-h AB u B } Y.rY., 



'(l + AY), 



(27) 
(28) 



where v, w and A are equal to the gauge invariant quan- 
tities of GS and GMG, again, as shown in Appendices lAl 
and[B] The vectors u A and n A indicate background unit 
vectors in the time- like and space-like radial directions, 
respectively. They are given by 



u A 



(1.0) 



and 



(29) 



The Hab m ([27|) correspond to the linear perturbations 
to gAB ! as shown in (|26p , and are included to ensure the 
normalisation u^u^ = — 1. 

The general form of axial perturbations to the metric, 
in RW gauge, are 



ds 2 =-dt 



(1 - nr 2 ) 
2k A Y b dx A dx b 



l —dr 2 + a\r 2 dU 2 (30) 
(31) 



where Ua is equal to the gauge invariant perturbation of 
GS, from Appendix |A1 Following GS, we also define for 
later convenience the new variable 



n 



AB 



kA 

a l r2 / W 



(32) 



1 We have changed notation from GMG to avoid potential confu- 
sion with notation that is often used in cosmology. In particular 
we have replaced: k ip, \ ip f , r\ i— » rj, a h-> v, (3 i— > 

v, w A, 7 h- > w. 
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where the fundamental anti-symmetric tensor eab — 
haUb — u A n b and a pipe denotes the covariant derivative 
on M 2 . The only axial perturbation that can then occur 
to the matter content is a perturbation to the matter four 
velocity, such that 



Scalars 



u M = (uA,vY a ), 



(33) 



where u A is defined as in the polar case, and v equals one 
of the gauge invariant variables of GMG, in Appendix [Bl 



1. 'Scalar-Vector-Tensor' variables 

The GI variables defined above give a very concise 
set of governing equations (see below). However, in the 
FLRW limit we shall see that they mix up the normal 
SVT modes in a complicated way - for example, the 
variable ip contains all three types of perturbations. By 
taking combinations of the variables defined in the per- 
turbed metric we may form variables which do reduce to 
scalars, vectors, or tensors in the FLRW limit, and so 
may be used to identify generalised SVT modes. Such 
variables are: 
Tensors 



polar: X = X. 



(34) 



axial: f = II" + 6WII' + ( 8W 2 — ^ +1 ) + 2 ] n 



-16tt 



(pv)> 



(35) 



Vectors 



polar: £ 



3a± 
2W 



V a i 



axial: v — v. 



(36) 
(37) 



\" + 2W\' - 



+ rWi' + r ( 3W 



2 9 



(38) 



where 



A = 8vrpa ± [H^A - 3u] . (39) 
We shall justify these variables in Sec. IIV CI 

III. MASTER EQUATIONS 

The equations of motion governing the dynamics may 
be reduced to a coupled system of evolution equations 
and constraints. For the polar sector the equations come 
in the form of three coupled PDEs for x, tp and with 
all other non-trivial gauge invariant variables determined 
from the solution of this system. For the axial modes, 
the dynamics are determined by a much simpler system 
of equations for the variables II and v. 

A. Polar perturbation equations for I > 2 



Substituting the polar perturbed metric tensor 
and the perturbed matter quantities, (j27|) and (j28|) . into 
the field equations results in a system of three coupled 
master evolution equations for the three variables Xi f 
and The remaining variables associated with the fluid 
can then be determined directly from the solution to this 
system. 

The three evolution equations are 



-X + x " -iH,X'2W x ' + 



(40) 



2(F|| - H±y - 2 F|j - 2(F|| - H±)W ? + i(H {l - H±)ip - 2 



8np- ^3- - 2H ± (H {] - H ± ) 



and 



(p + 4H ± p - 2 — W 2 ) ip = -H x x + W X ' 



2W 2 



£(£+1) + 2 



2a\r 2 



X +2W(H ll -H 1 _)<;, (41) 



and 



? + 2if,K = -x' 



(42) 
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together with the constraint 



f] = 0. 

These four equations are the LTB version of Equations (GMG87), (GMG88), (GMG89) and (GMG82) from @. 
The three constraint equations can then be written as 



(43) 



8itpw 
8irpA 



- (H\\ - 2tf x y -Wx + H ±X ' 



2a 2 r 2 

-<p" - 2W<p' + (H\\ + 2H x )<p + W X ' + H±x 
(l-W + 2) 



£(£+!) + 2 , , 

' f Hi + 2H ± Hn -W 2 - 4np 



1 2Hl + AH ll H 1 _ ~ 8np 



a^r 2 



2a 2 r 2 



-X + 2H ± q' + 2(H«+H X )W<; 



X <;' 
8npv = <p + - +H\\(x + <p) + -z 



(44) 

(X + V) (45) 
(46) 



These are the LTB versions of the equations (GMG93), (GMG94) and (GMG95) from |||. It is also useful to consider 
the evolution equations that result from differentiating (I44|) . (|45l) and (|46|) . These take the particularly simple form 



X . 9 

v = , 

2 2' 



w 
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A 



#11 

Huw h, 

2 11 2 ' 

X 3y £{£ + !) _p^_ 
2 2 a\r 2 V p 



2W(w 



(47) 
(48) 
(49) 



and are the LTB versions of (GMG96), (GMG97) and 
(GMG98). There are then six coupled equations for the 
six variables x, ip, w, A and v. 

Equations ([4^]) - (j4^)) are clearly more complicated than 
their FLRW counterparts, but the fact that they can be 
written as concisely as they are above is quite remarkable. 

For very large angle fluctuations, with £ = or 1, the 
field equations no longer give 77 — 0, which has been 
used to simplify the equations in this section. Instead, 
there are additional gauge freedoms that can be used to 
simplify (|4"TJ)) - (|49p . These are discussed in Appendix [Cl 



B. Axial perturbation equations 

The axial perturbation equations take a simpler form 
than their polar counterparts. Because we are consider- 
ing dust dominated cosmologies, the field equations give 
us v — 0, and so v — v(r) and must be set by initial 
conditions. For £ > 2 the metric perturbation LI can be 
shown to obey the wave equation 



-ii + n" - (6H± + #n )ri + 6\vn' 



(50) 



16tt 



(pv)> 



which is the LTB version of (GMG69). This is the only 
dynamical equation that needs solving, for the single vari- 
able n. Once this equation is solved it gives the axial 



metric perturbations, k>A, via 

{t-T)(£ + 2)k A = \<a-K P a\r 2 vu A - <t AB (a\r^ B . (51) 

Clearly £ = 1 is a special case. For large angle perturba- 
tions of this type we have (a^r 4 !!) = and {a^r^Ii) 1 = 
— Wira^r 2 pv. The solution for n is then 



n 



v(Mr 3 ), r 
a\r 4 J VI - nr 2 



dr. 



(52) 



The metric perturbations k A must then be obtained from 
inverting Eq. J32|) . This will involve an extra degree 
of freedom in k A that means it cannot be determined 
uniquely. 

There are no axial perturbations with £ = 0. 



C. Structure of the equations and solutions 

It can be seen from the master equations, (14T)|) - (14"9"1) . 
that x an d n contain gravitational wave degrees of 
freedom, as their leading derivatives are of the form 
— ( ") + ( )" (i- e - a wave equation with a characteris- 
tic speed of unity). It can also be seen that ? and v look 
rather like vectors in FLRW cosmology, while ip appears 
to govern density perturbations like the Newtonian po- 
tential. Unfortunately, this roughly analogous behaviour 
does not hold completely, as we will show below when we 
consider the FLRW limit. What is clear, however, is that 
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there exist complicated couplings between the variables 
in this system. We can therefore no longer expect gravi- 
tational waves to decouple from density perturbations, as 
the varying curvature of the background serves to couple 
these perturbations intimately. 

The solutions to the perturbation equations in sections 
IIII Al and IIII Bl will, in general, need to be found numer- 
ically. However, if we consider i > 2 and neglect the 
contribution of the gravitational waves described by Xi 
we can make progress. With 77 = x = 0, Equation 
can be integrated to find 



C oc 



1 



Clearly, <; — > as the Universe expands, and au — > 00. 
Assuming \ to be negligible can then be seen to lead to 
a negligible <j in the late Universe. With \ = ? = 0, 
Equation (HT]) becomes 



2k 

<f + 4Hj_ip T ip = 0. 



(53) 



This equation can be solved parametrically, together with 
the analogue of the Friedmann equation, ([3]). There are 
three solutions, depending on the sign of k. 
For k < we find the solution 



aj_ = 



M(l - cosh 29) 
2k 



(54) 



cosh (P) 

tp = . \ 5 [c 1 +c 2 (sinh29-6e (55) 



sinh 6 



-4tanh6)] 



t-tr 



M(sinh29 - 29) 

2(-k) 3 / 2 : 



(56) 



where C\ = c x (r), c 2 = c 2 (r) and t = t (r) are free 
functions that must be specified as part of the initial 
conditions. The function to(r) is the 'bang time'. 
For k > we find 



a_L 



M(l-cos29) 
2k 



■_- = C ° S ® [ci + c 2 (sin 29 — 69 



t-to = 



sin 9 



M(29-sin29) 



+4 tan 9)] 



2(k) 3 / 2 



(57) 
(58) 

(59) 



Again, ci , c 2 and to are free functions of r. 

Finally, for k = 0, a parametric solution is not re- 
quired, and we find 



(18M)V3( t _ t ) 2 / 3 



V = 



ci 



c 2 . 



{t-tof/z 

Once more, with ci, c 2 and to as free functions of r 



(60) 
(61) 



In the c / cases, 9 is a monotonically increasing 
function of t, and in all three cases c\ corresponds to a 
decaying mode, and c 2 to a growing mode. In the cases 
with k < the 'growing mode' is itself either a constant 
(for k = 0) or a decreasing function of t (for k < 0), and 
in those cases we mean growing with respect to the other 
mode, which is decaying faster. 

For k — the power-law forms of a±(t) and if mean 
that the asymptotic form of both the growing and decay- 
ing modes is clear. It can also be seen for k ^ that 
both modes, together with a±, behave as if k = in the 
limit t — ► to (as expected, due to their dynamical evolu- 
tion being dominated by p, and not k, at early times). 
In the case of n < the growing mode can be seen to 
decay as ~ l/(f — to) when t — > 00, and the scale fac- 
tor, a±, behaves like an open Milne universe. For k > 
the scale factor initially grows, reaches a maximum of ex- 
pansion, and then collapses to zero in finite time, t. As 
this future singularity is approached, both growing and 
decaying modes of tp diverge to infinity. 

These solutions may serve as the basis for calculating 
the matter power spectrum. However, they will only re- 
main so as long as the evolution equation for \ remains 
approximately satisfied. 



IV. THE FLRW LIMIT 

From the perturbed LTB equations above it is not 
clear what each gauge invariant variable corresponds to 
in terms of the more familiar Bardeen potentials, and so 
forth, of standard cosmology. Here we derive the stan- 
dard FLRW perturbation equations in the GMG formal- 
ism, and re-express the GS and GMG variables in terms 
of FLRW variables. 

In the FLRW limit k — > constant, and we have that a± 
and <z || — > a(t). The master equations in the polar sector 
(gni)-(I121) then become 



dl + 2Ud T - V 2 



A(l-Kr 2 ) d 2 



x = o, 



and 



[d 2 + md T - 2k] <f 

_ ndT+ (l^ dr+ W + l)-2 + ^) 



2r 2 



and 



[d T + 2H]s= v 7 ! ~ Kr 2 d rX , 
while in the axial sector we have 



d 2 T + 6Hd T - v 2 - 4(1 Kr } f d r + ^ ) + m 2 



3 

2r 



n 



Vl - Kr 2 _ 
lbirp ^ o r v, 



(62) 
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and 



d T v = 0, 



(63) 



Let us now consider the perturbed FLRW line-element 
in the longitudinal, or conformal Newtonian, gauge: 



where dr = dt/a is conformal time, and TL = a. T /a is the 
conformal Hubble rate which obeys the Friedmann and 
Raychaudhuri equations 



ds 2 



-a 2 (I + 2$)dr 2 - 2a 2 V l drdx l 



(65) 



[(1 - 2¥)7* 



hij\dx l dx\ 



H 2 = ^-k 



and d T H= - ~(H 2 + k) . (64) 



Throughout this section V 2 will always refer to the Lapla- 
cian acting on a 3-scalar, so that 



,> , (2-3^r 2 ) g _ £(£_ 



1) 



V 2 = (1 - nr 2 )d; 



It is evident from the equations above that x is a grav- 
itational wave, as its characteristics are null. However, 
X can also be seen to act as a source for (p even though 
the homogeneous part of the evolution equation looks 
like very similar to the Bardeen equation. Similarly, s 
looks like it almost obeys the vector decay law, but is 
also coupled to gravitational waves through \- The fact 
that these variables do not decouple in the FLRW limit 
means that their interpretation in terms of FLRW gauge 
invariants will not be straightforward. 



where a = a(r) is the scale factor, and 7y = diag[(l — 
rer 2 ) -1 , r 2r y a b] is the background spatial metric in spheri- 
cal coordinates of curvature n. V, is the covariant deriva- 
tive with respect to jij. The metric (|65[) is split in the 
standard way into two 3-scalars ($, a 3- vector (Vi) 
and a 3-tensor (hij), where the '3' is not usually empha- 
sised. All of these quantities are gauge invariant variables 
(see Appendix |D|) . The vector Vi is divergence-free, and 
the tensor is divergence and trace-free. Note that the 
coordinates used in the metric in the longitudinal gauge, 
given by Eq. ([55]) , and the metric in the RW gauge, Eqs. 
(f26l) and (1301). are not the same. 



In Appendix iDl we show that, in the FLRW limit, the 
LTB gauge invariants (<p, Xi Vj EE and v) can be written 
in terms of the usual FLRW invariants ($, 'J, Vi and hij) 
in the following way: 



J 



Polar: 



V = 



-29 - 2HV - 2 (1 Kr \ r + ^h< T > 



^ + (l_^)^ + ^ + l)-4(l- K . 2 ) 



X = (1 — nr 2 )h rr + 2 



-(1 - nr 2 )d r + 



2\n2 2nr 2 ) 
(1 — nr )o r o r — 



2r 2 



y/l - kt 2 |y r - 8 r V + d T h r - (d r - ^ 9 r /i rrF) I 



K - —h (T) 

+ 1) - 8 + 4nr 
2^2 



2 1 



(TF) 



$ - * - (d T + 2H) V + -(1 - nr 2 )h rr + [-(1 - nr 2 )d r + ttr] h r 



-d 2 + (1 - Kv 2 )d 2 r - 2-Hd T - (2 Kr2) a r + 4 



(66) 
(67) 

(68) 
(69) 



Axial: 



n 

Wirpav 



s/T 



d r ] V + d T h r 



-V 2 + ^ -d r - 4k 



r 



(70) 
(71) 



where, in order to compare with the GMG formalism, 
we have split the 3-vector Vi and the 3-tensor into 
their radial and angular parts, and then these into their 
polar and axial components, as described in Appendix 
[Dl This makes explicit the mixing of SVT modes in the 



LTB gauge invariants. We have substituted Eq. (|D10[) to 
remove the tensor contribution to v in Eq. (|7ip . 



A. Polar Perturbations 

As we are dealing with FLRW perturbations we may 
separate the scalar, vector and tensor parts of these equa- 
tions, as we know they evolve independently. 
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1 . Scalars 

From 77 = we have $ = $ , as is usual in a dust 
dominated FLRW cosmology. From the (/j-equation we 
then find 

<9 2 * + 3Hd T y ~ 2k* = 0, (72) 

which is the usual Bardeen equation for the Newtonian 
potential. The <;- and ^-equations contain no scalar 
modes. We then note that the scalar part of the gauge 
invariant matter perturbations can be written 

47ra 2 pA = V 2 * - md T ty - 3(H 2 - (73) 
A-Kapv = -{d T +H)V, (74) 
4wa?pw = -y/l - nr 2 {d T + H)d r V. (75) 

It can be seen from ([75)1 that the scalar part of the energy 
density perturbation is just the gauge invariant density 
fluctuation Sp (G " = Sp + d T p(B — d T E) However, 
we will see below that we cannot simply identify A = 
5p <GI, /p, as A has non-zero vector and tensor parts. 

2. Vectors 

From r\ = 0, and the evolution equation for <;, we find 
that 



of the FLRW perturbed metric, as the x equation con- 
tains fourth derivatives of h <TF) . In fact, as all four GMG 
variables contain 3-tensors, all four field equations con- 
tain tensor modes, and showing that they reduce to the 
usual FLRW wave equation for tensors is non-trivial. The 
four equations form a set of coupled PDEs for the vari- 
ables h rr , h r , h <T) , h <TF) , the simplest of which may be read 
off from 77 = above. Recall, however, that there is only 
one degree of freedom in the polar part of hij, with the 
trace and divergence free conditions removing the other 
three. If we choose this degree of freedom to be h (T> then 
it can be shown from the GMG equations for 77, tp and f, 
when combined with the conditions given by Eqs. (|D7[) . 
(|D8|) and (|D9|) . reduce to the single master equation 

{d 2 + 2Hd T -V 2 }h m = 0, (80) 

which is just the trace of the usual FLRW wave equation. 
The x equation then follows automatically. 

The tensor parts of the variables associated with the 
fluid can then be written 

A = hidrh™, (81) 

v = ^ad T h <TF) , (82) 

w = ly/i- K r 2 (d T h r + ^d T h (T A . (83) 



d T V = -2HV and d T V r = -2HV r . 



(76) 



Axial Perturbations 



Thus, these two equations give the usual a~ 2 law for 
vector modes. The vector part of the evolution equation 
for ip automatically follows from these two, while the x~ 
equation contains no vectors. The vector parts of the 
gauge invariant matter perturbations are then given by 



A 

I6irapv 
I67ra 2 p w 



3HV (77) 
[(1 - nr 2 )d r - nr] V r (78) 
- [(1 - Kr 2 )d 2 r - nrd r - 2{m 2 + k)] V, 
'£(£+1 



£(£+!) 



+ 



3{H 2 + k) 



0. 



It is a curiosity that the vector part of A = 3HV is non- 
zero, so that while A looks like it is just a gauge invariant 
density perturbation, it actually contains vectors. 



3. Tensors 

The GMG formalism is rather less tidy when trying 
to describe 3-tensor modes. The equation which looks 
like it should be describing gravitational waves is that 
for x- But while this equation consists only of tensor 
modes, it is in a very ugly form when evaluated in terms 



There are no scalar modes in the axial sector. 

For vector modes it can be seen from d T v = 0, and 
the radial equation of l[5l|). that the usual vector decay 
law, d T V = —2HV, is obeyed. The Il-equation is then 
automatically satisfied. 

Combining the radial equation of (|?T|) and Eq. (jDlOl) . 
we find the wave equation 



[d 2 + 2Ud T - V 2 + 2nrd r + 3n]h r = 0, 



(84) 



V r (79) where the extra curvature terms are due to V acting 
here on a scalar that is actually part of a 3-tensor. The 
equation for LI then follows. 



C. SVT Variables 

We have found that what appear as gauge invariant 
fluid perturbations in the GMG formalism - namely A, 
v and w — are not exclusively fluid modes, as they are 
usually understood in FLRW cosmology, as they can be 
excited by tensors (FLRW models need a tensor part in 
the anisotropic stress for fluid perturbations to couple to 
tensor modes, which we have not included here). It is 
also the case that the GMG metric perturbations do not 
correspond in a straightforward way to the SVT decom- 
position we are familiar with from FLRW cosmology. For 
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example, it is clear that x represents gravitational waves, 
and in the FLRW limit it is a pure tensor mode: but what 
about scalars and vectors? Can we identify combinations 
of gauge invariant GMG variables which will represent 
purely scalar or vector modes in the FLRW limit, and so 
be useful physical variables in LTB cosmology? 

In order to find these, let us first note that the combi- 
nation 



A = 8na 2 H~ 1 p 



A - 3—w 

a 



(85) 



does not contain any tensor modes in the FLRW limit (as 
can be verified from Eqs. (fSTj) and ([52])). Also, note that 
we can construct from ^ and \ a quantity that contains 
only vector modes: 



3r 



2^fT^ Kr 
V 2 + :!/v 



[V 2 +3K]r-\+-^d r (r 2 d TX ) 



~-d r 



(V r ~ 8 r V) . 



(86) 



From A and £ we can now construct a variable that is a 
function of scalar modes only: 

3nr 2 ' 



C 



a~ 2 V 2 A 



(1 - nr z )d r + 



a 2 H 



3k 



V 2 *. 



(87) 



This is related to the gauge invariant density fluctuation 
8p (a,> and the curvature perturbation, 1Z = W — H(H& + 

d T ^)/{d T n-n 2 -K) pa, by 



2K 

8np 



= V 2 



5p<' 



3ft -3* 



L P 



For the axial perturbations we already see that v is a 
pure vector mode, but that II is an awkward mixture of 
vectors and tensors. Defining the variable 



T = 



- a- 2 



^ 2 + i(l-nr 2 ) d 2(3 - 4kH) 



n 



16np 



-d r v 



vT 



v 2 



2nrd r 



(89) 
(90) 



we are able to define a pure tensor mode using the GMG 
gauge invariant variables. 

Using these new variables, we can now generalise to 
the inhomogeneous case. The resulting functions will 
have the useful feature that they reduce to pure scalar 
or vector modes in the FLRW limit. An example set 
of variables is given by Eqs. (j3~i|) - (f3"9")) . in section Hi C II 
These reduce to the quantities above in the FLRW limit. 
The variables £> v, x an d T, are then useful gauge 
invariant variables in the LTB perturbation equations as 
they encode generalised scalar, vector and tensor pertur- 
bations, respectively, as they are normally referred to in 
cosmology. 



V. DISCUSSION 



We have presented here a full system of master equa- 
tions that represent the general perturbations to LTB 
space-times, in terms of the gauge invariant variables 
ip,<;,X, II, v (and rj, strictly speaking). This formalism 
can now be used to investigate, in a fully consistent fash- 
ion, the growth of linear structure in LTB models. As 
such, predictions for phenomena such as baryon acoustic 
oscillations, the integrated Sachs-Wolfe effect, and the 
observed matter power spectrum can now be made, and 
used to compare with astrophysical data. This will allow 
us to establish what deviations should be expected from 
the usual FLRW predictions, and whether or not such 
deviations can be observed in our Universe. 

However, while <p,c, Xi-H, v and r\ appear completely 
naturally, and produce an elegant system of equations in 
the form of an initial value formulation, they do not re- 
duce to anything intuitive in the FLRW limit. Instead, 
they are a cumbersome mixture of scalar, vector and ten- 
sor perturbations - a situation made worse still by the 
additional couplings that exist in the evolution equations 
in the more general, inhomogeneous case. 

To address these problems we have defined a new set 
of variables in the polar sector, £, £ and x, that can 
be thought of as generalised scalar, vector and tensor 
perturbations, respectively. Similarly, we propose new 
variables in the axial sector, v and T, that correspond 
to generalised vector and tensor modes. In the FLRW 
limit these new variables become pure scalars, vectors 
and tensors. We expect these new variables to prove use- 
ful in future studies where they can be used to set the 
initial conditions for the evolution of perturbations. For 
example, in an LTB model with homogeneous bang time, 
to(r) ^constant, the surfaces of constant t will be almost 
homogeneous at early times, and will then be well ap- 
proximated by an FLRW description. The new variables 
can then be used to match the FLRW initial conditions 
onto the LTB evolution equations in a straightforward 
way. This will be considered further in future studies. 



APPENDIX A: THE GERLACH-SENGUPTA 
FORMALISM 



The general form of perturbations to the metric and 
stress-energy tensors can be written, in terms of the 
harmonic functions outlined in Section HU as — > 
g^u + Ag^ and T^„ -> T^ v + AT^ V , where [TO] 



Aguu = 



AT,,. 



_ hf i3l Y a 
hT^Y a h Y ab 

o A^ xial y Q 

Atf ml Y a At^Y ab 



(Al) 
(A2) 
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for axial perturbations, and 



Ag^ = 



AT,, 



lAB 
3olari 



Y 



, polar v 
IL A a 



h p A olaT Ya a\r 2 {KY lah + GY.. ab ) 



At AB Y 
At p / ai Y a At&Y a 



(A3) 



a .polar 
A 



Y, 



i\r 2 At^Y lab 



(A4) 



for polar perturbations. The new variables hjy, h, 

AtT*\ Ai(D, hAB) Kj Gj MaBi At Polar ; At(2) 

and At^ 3 ) are all functions of the two coordinates x A . 

From these quantities we can then construct the gauge 
invariant variable [Till 



k A 



u axial 

n A 



h\ A + 2hv A , 



for axial perturbations, and 
kAB = h AB - Pa\b - Pb\a 



ip = K - 2v A p A 



t A B = AtAB 

t A = AtP° la 



t A B\CP 
~ t AC p C 



tACP C \B 



tBCP C \A 



The quantities At^ 1 



(A5) 

(A6) 
(A7) 
(A8) 
(A9) 

At(D, 



for scalar perturbations. _li±c quanuuco L -* b A 
At^ and At^ are already gauge invariant for a dust 
dominated universe. Here a pipe indicates a covariant 
derivative with respect to the metric on M 2 , and we have 
defined pa = h A ° Ur — a 2 L r 2 G\ A /2, along with GS, for con- 
cision. All of these expressions are constructed so as to be 
invariant under infinitesimal coordinate transformations 
of the form x^ — > x^ + £ M . 

There exists here a useful gauge in which h = ft.^° lal = 
G = 0. This is the Regge- Wheeler gauge [15] . where the 
remaining perturbation variables are all equal to gauge 
invariant quantities, as can be seen from the expressions 
above. 



APPENDIX B: THE 
GUNDLACH-MARTIN-GARCIA FORMALISM 

The background quantities used by GMG can be writ- 
ten in the LTB case, using our notation, as 



GMG then write their perturbed fluid four velocity as 
u n u [t + Au M , where 



Au M = (0, v Y a ) 
for axial perturbations, and 



Aw,, 



wn A + -h A BU B )Y,vY a 



(B6) 



(B7) 



for scalar perturbations. Here ha = ~~^abu b is a unit 
space-like vector, and v, v and w are all functions of 
x A . Density perturbations are then parametrised by p — > 
p+AYp, and the following gauge invariant quantities can 
be constructed: 



V — V — p U B 



w = w — n A UA\BP B 



A ee A-p A {\np\ A 



-n A u B 



(B8) 

(pa\b-Pb\a) (B9) 
(BIO) 



with v already gauge invariant. Again, the Regge- 
Wheeler gauge can be seen to have a special significance. 
The gauge invariants above can be related to those in 
Appendix [A] via 



L A = vu A 
t A = vpu A 



tAB = P 



w(uAriB + tiaUb) + Auaub 



+ -ji^ACUB + U A k B c) 



(Bll) 
(B12) 

(B13) 
(B14) 



together with L = < (2) = i (3) = 0. GMG then continue 
to decompose kAB into the 'fluid frame' via 

kAB = T}{n A n B - uaub) + 4>{n A n B + u A u B ) (B15) 

+ s(u A n B + n A u B ), (B16) 



u = 


H± 


(Bl) 


V = 


H± 
W 


(B2) 




M 


(B3) 


m = 


~2 


I 1 = 


H \\ 


(B4) 


v = 


0. 


(B5) 



and to define the new variable 



if + T]. 



(B17) 



In these variables the perturbation equations take a par- 
ticularly simple form [8[. 
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APPENDIX C: POLAR i = 0, 1 PERTURBATIONS 

For polar perturbations with I = 0, 1 it is no longer the case that the field equations give rj = 0. The general system 
of equations given in Sec IIII Al now reads: 
For l>\: 



= -2 



(CI) 



3M 

8np- — -2H ± (H ll -H ± ) 



(x + p) + {£ y 2 +2) X + 3H\\x + 4(H|| - ff ± )y> + 2W X ' 



H' - 2(H\\ - H ± )W 



— Cp = 4H±ip 



<;-2(H l{ - H±)f] - 6Wr)' - 
2 £(£ + !) + 2 



£(£ + !) + 8 



8H l] H ± + AH'i ~ 16W 2 - 32irp 



2W Z 



2a\r 2 



( X - 277) - 2W(H« - Hx)s 



1 



2 2 



+i?±(x - 2?)) - WOC - 277') - 2W 2 r] - 2 



i = 2H^ + x ' + 2Wt}-2if. 



X . <P 
v = — 77, 

2 2 '' 



— W ip, 



(C2) 

(C3) 
(C4) 
(C5) 



and for i > 0: 



87rpw = (0)' - + F_lx' - (-H"|| - 2H x )<p' - 2H x r)' 



*V±V+± + rrl 1 2// // sr-' ■ -i^ 



8?rpA = 



IU + 1) 2 

V 22 + 2Hl + 4H\\Hx - 8np 



(x + p) + 2Hx<; 



, (t-W + 2)_ 



A+ r + 2J =^i^ v -r + 2J7"2-T- 2iy r + 2J 



2a^r 2 X 

+2(ff|| + H±)Ws + H ± x - p" + {H\\ + 2H ± )ip + W X ' - 2Wp' - 2H ± (H ± + 2H\\)n, 
w = — - H\\w - — <r - W-q, 



<:, (C6) 
(C7) 

(C8) 
(C9) 



1. 1=1 



For I = 1 there is an additional gauge freedom which can be used to eliminate one of the perturbation variables. 
The obvious choice may be to use this to set 77 = 0, but GMG showed this to lead to some ambiguity. We therefore 
follow them in using the less ambiguous choice p> = 0. The evolution and conservations equations above can then be 
combined to give 

3W 2 + H'i J 77 - [W 2 + Hi - 4ttp] x - 2H ± W<;, (C10) 



Wrf - 


HxV 


w x '- 


Hxx 


wV- 





2 



2H± ~ 8irp 



Wq'-Hx<i = 8Trp(w + 2Wv)+2HnW(T)-x) + 4HxWT) 



X -2(H l{ +Hx)W<;-2Hlr) 
2 



Hi - W 2 



Airp 



(Cll) 
(C12) 
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These are the LTB equivalents of (GMGA12), (GMGA13) and (GMGA14) from @. These three equations can be 
solved, along with the three conservation equations (|47[) - (|49p . with (p = 0, for the six remaining perturbation variables, 
77, Xj v , w and A. 



= 



Polar perturbations with I = are spherical. As such one could conceive of absorbing such perturbations into the 
background parameters, which themselves describe the most general spherically symmetric dust dominated solution 
of Einstein's equations. Despite this being the case, it still seems worth considering I = perturbations in a spherical 
background, as one can then separate them out more easily from a smooth background. 

In this case we again have a non-zero 77, but now have sufficient gauge freedom to set tp = and 



2H 1 _W 
Hl + W< 



:(V - X), 



as done by GMG. The constraint equations then give 
Wr)' — H±fj = Airp I \ 



(W 2 



I6irp 



{W 2 - HI) 



A {W 2 + Hi) 
2 w + top 2 , A, 



(W 2 - HI) 



and 



W X '-H ±X = 



AH ± W 
(Hi + W 2 ) 



H±W + 4-717? 
H 2 



H L W 



(W 2 - H\ 



-A 



%np- 



(W 2 - Hi 
1 



(x-v) + 16?r P 



(W 2 - Hi) 



(C13) 



(C14) 



(C15) 



2H 2 ± 



(W 2 - Hi) 



These are the LTB versions of (GMGA18) and (GMGA19). The two remaining conservation equations are then 



- H\\ + 4irp 

+(v-x) 



H t 



(W 2 - H 2 ± [ 
W{W 2 - Hi 



W 



2a\r 2 {W 2 - Hi) 



(C16) 



H\\H±W 



2(W 2 + Hl) (W 2 



Hi) 



4irp 



WH 2 



and 



A + w 



iwpHj 



(W 2 - H 2 ± ) 



-(y-x) 



A 



H 



(W 2 - H\ 



2a\r 2 



4irp r\ — 



{Hi - W A ) 



2W - Airp 



W 



(W 2 - Hi) 



(C17) 



H±W 



(W 2 + Hj 



H ± +H 



H X (W 2 - Hi 
2(W 2 +H 2 ± ) 



Airp 



H ± W 2 



(W 4 - Hj) 



These equations are the LTB versions of (GMGA15) and (GMGA16). We now have four equations to solve for the 
four remaining variables 77, x, w and A. 



APPENDIX D: MATCHING LTB AND FLRW 
GAUGE INVARIANTS 



Here we derive the standard FLRW perturbation equa- 
tions in the GMG formalism. Unfortunately, the RW 
gauge is not well adapted to our usual description of 
FLRW perturbations. We must therefore write both the 
GMG variables and the perturbed FLRW metric in a 
general gauge. 

In general coordinates an arbitrary perturbation of 



FLRW can be written 

ds 2 = -a 2 (l + 2(j))dT 2 + 2a 2 (ViB - S^drdx* (Dl) 
+ a 2 [(l - 2ip)jij + 2V l V j E + 2ViFj + h ij ]dx i dx j , 

where a = a(r), 7^ is the spatial metric, and Vj is the 
covariant derivative with respect to 7y . The vectors here 
are divergence free, ViS* 1 = 0, and the tensors are diver- 
gence and trace-free. Gauge invariant metric perturba- 
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tions are given by 

$ = <f> + H(B-d T E) + d T (B-d T E), 
* = ip — TL(B — d T E), 

V, = Si + d r Fi, (D2) 

and the perturbations hij are already gauge- invariant. 

To compare with the GMG formalism, we expand this 
1+3 split into a 1+1+2 split in spherical coordinates, 
peeling off the radial parts of each variable. A 3-vector, 
such as Si, then splits into a scalar, S r , and a 2- vector, 
S a . These can then be decomposed into spherical har- 
monics as 

S r = J2 S r m)Yiem) = S r Y , 
£rn 

g _ g(fm)y-(£m) ^_ g(^ m ) e b Y^ m ^ 

= SY a + SY a , (D3) 

where S and S are the polar and axial parts of S, re- 
spectively The divergence-free property of the 3-vector, 
ViS* 1 = 0, then gives us that 

(1 - nr 2 )d r S r + {2 ~ 3 r Kr2) S r - fc^S = 0. (D4) 



This reduces the degrees of freedom in the 3-vector to 
two: one polar and one axial. 

For a 3-tensor, the split into radial and angular parts 
is a bit more messy. The tensor 





h ra 


h ra 


h a b 



(D5) 



splits into a 2-scalar h rr = h r e r m ^Y^ m ^ (spherical har- 
monic sum implied), a 2- vector h ra = h r lm ^Ya im ^ + 
^(£m)y^m) a remamm g p ar ^- cons isting of a 2-scalar 
(the trace) and polar and axial trace-free 2-tensors: 

h ab = h' T > lab Y + h' TF >Y ab + hY ab . (D6) 

This splits the 3-tcnsor hij into its polar (h rr , h r , 
and h( TF " > ) and axial (h r and h) parts. 

The trace- free, 7*- 7 'ft»j = 0, and divergence- free, 
V % hji = 0, conditions give: 



= (1 - nr 2 )h rr + \h (T) , (D7) 



= (1 - Kr 2 )d r h rr + (2 W) - Xh m - ^ % l) hr, (D8) 

= (1 - nr 2 )d r h r + (2 ~ 3Kr2) h r + \h (T > - V-W + V hW (D9) 

r r z 2r A 

= (1 - nr 2 )8 r h r + {2 -^h r - «-W + 2 k. (D10) 



Again, the degrees of freedom left after applying these relations is two: one per parity. 

We can now equate the perturbed FLRW with the perturbed LTB metric, in an arbitrary gauge, to find that the 
polar perturbations are related by 

a -2 h GMG = _ 2 ^ (Dn) 

a- 2 h r : r = - T -^— s il> + 2 (V - T ^-^] d r E + 2(d r - — ^— F r + h rr , (D12) 
1 — nr z \ 1 — kt z J \ 1 — Kr z J 

a~ 2 h^ G = -S r + d r B, (D13) 

a -2 h GMG = _ S + Bj ( D14 ) 

a -2 h GMG = 2^d r -^jE + F r +(d r - 2 ^jF + h r , (D15) 

r 2 G = 2E + 2F + h (TF> , (D16) 

K = -2yj + 2 X —^d r E + 2 X —^F T + \h m + l -^fi-h< TF \ (D17) 
r r r A 2r 

where everything has been decomposed into spherical harmonics, so that <j> — <p^ m '(x A ) etc. For the axial modes we 
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have 



a- 2 h° MG 


= -s, 




a~ 2 h GMO 


= (d r 


2 

r 


ar 2 h GUG 


= F + 





(D18) 
(D19) 

(D20) 



The gauge invariant GMG variables can now be calculated directly, and are given in Section IIV 
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